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Abstract 

Background: Prion diseases are characterized by the accumulation of the pathogenic PrP Sc protein, mainly in the 
brain and the lymphoreticular system. Although prions multiply/accumulate in the lymph nodes without any 
detectable pathology, transcriptional changes in this tissue may reflect biological processes that contribute to the 
molecular pathogenesis of prion diseases. Little is known about the molecular processes that occur in the 
lymphoreticular system in early and late stages of prion disease. We performed a microarray-based study to identify 
genes that are differentially expressed at different disease stages in the mesenteric lymph node of sheep naturally 
infected with scrapie. Oligo DNA microarrays were used to identify gene-expression profiles in the early/middle 
(preclinical) and late (clinical) stages of the disease. 

Results: In the clinical stage of the disease, we detected 105 genes that were differentially expressed (>2-fold 
change in expression). Of these, 43 were upregulated and 62 downregulated as compared with age-matched 
negative controls. Fewer genes (50) were differentially expressed in the preclinical stage of the disease. Gene 
Ontology enrichment analysis revealed that the differentially expressed genes were largely associated with the 
following terms: glycoprotein, extracellular region, disulfide bond, cell cycle and extracellular matrix. Moreover, some 
of the annotated genes could be grouped into 3 specific signaling pathways: focal adhesion, PPAR signaling and 
ECM-receptor interaction. We discuss the relationship between the observed gene expression profiles and PrP Sc 
deposition and the potential involvement in the pathogenesis of scrapie of 7 specific differentially expressed genes 
whose expression levels were confirmed by real time-PCR. 

Conclusions: The present findings identify new genes that may be involved in the pathogenesis of natural scrapie 
infection in the lymphoreticular system, and confirm previous reports describing scrapie-induced alterations in the 
expression of genes involved in protein misfolding, angiogenesis and the oxidative stress response. Further studies 
will be necessary to determine the role of these genes in prion replication, dissemination and in the response of 
the organism to this disease. 
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Background 

Transmissible spongiform encephalopathies (TSEs) or 
prion diseases are a group of fatal neurodegenerative dis- 
orders of animals and humans. TSEs generally involve 
long incubation periods, which in humans can span sev- 
eral decades [1]. These diseases are characterized by the 
accumulation, mainly in nervous and lymphoreticular 



* Correspondence: rbolea@unizar.es 

Centra de Investigacion en Encefalopatias y Enfermedades Transmisibles 
Emergentes, Facultad de Veterinaria, Universidad de Zaragoza, Zaragoza, 
Spain 

Full list of author information is available at the end of the article 

(3 BioMed Central 



tissues, of an abnormal isoform (PrP c ) of a normal 
host-coded cell surface glycoprotein (PrP c ). Pathogenesis 
typically involves a triad of histological events; variol- 
ation, neuronal loss and astrocytosis. PrP Sc accumulation 
is currently considered the only reliable biochemical 
marker for this group of diseases [2] . 

Scrapie is a TSE that naturally affects sheep and goats. 
Ovine scrapie is a useful model in which to study the 
pathogenic mechanisms of the variant Creutzfeldt Jakob 
disease (vCJD), which affects humans [3]. Like vCJD, 
classical scrapie is associated with widespread tissue in- 
fectivity, mainly in the central nervous system (CNS), 
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peripheral nervous system (PNS), and lymphoreticular 
system (LRS) [4]. 

Mesenteric lymph nodes are one of the first tissues of 
the LRS in which PrP Sc accumulates during the preclin- 
ical disease stage in naturally infected sheep [5,6], al- 
though this process strongly depends on the animals 
PRNP genotype. Subsequently, other LRS tissues are rap- 
idly and near simultaneously exposed to infectivity, pre- 
sumably by circulating prions. In the clinical stages of 
the disease and in all susceptible PRNP genotypes except 
VRQ/ARR, infectivity and PrP Sc accumulation are gener- 
ally observed throughout all lymphoid tissues, except for 
the thymus [7]. Follicular dendritic cells (FDCs) may be 
necessary for prion propagation within the LRS [8], and 
macrophages and tingible body macrophages (TBMs) of 
the LRS have been identified as reservoirs of the infec- 
tious agents of TSEs [4,9]. 

The molecular mechanisms that underlie the uptake of 
the infectious prion agent and the progression of the dis- 
ease remain largely unknown. Given the complexity and 
multifactorial nature of the spread and accumulation of 
the infectious agent, we and others have used gene ex- 
pression analysis platforms to identify signaling pathways 
that are altered during PrP Sc accumulation and subse- 
quent neurodegeneration. This approach facilitates the 
identification of potential biomarkers and drug targets in 
natural ovine scrapie [10,11] and in experimental murine 
scrapie models [12-14]. 

Few gene expression-profiling studies have investigated 
changes in the lymphoreticular system in sheep with nat- 
ural scrapie; most have been performed using murine 
models [15,16] or ileal Peyers patches from orally inocu- 
lated lambs [17]. In the present study, we used an ovine 
microarray technology to identify transcriptional changes 
in the mesenteric lymph node in both clinical and pre- 
clinical disease stages in sheep that were naturally in- 
fected with scrapie. This is the first transcriptome-wide 
expression profiling study of the lymphoreticular system 
in sheep with natural scrapie. 

Methods 

Ethics statement 

This study was approved by the Ethics Committee for 
Animal Experiments of the University of Zaragoza (Permit 
Number: PI02/08) and was carried out in strict accord- 
ance with the recommendations for the care and use of 
experimental animals and in agreement with national law 
(R.D. 1201/2005). 

Animals, necropsy, tissue collection and RNA isolation 

All details of the animals used and the necropsy and 
RNA isolation procedures have been previously reported 
[10,11]. Briefly, we used a 17 female ARQ/ARQ Rasa 
Aragonesa sheep (aged 1-6 years). Eleven animals were 



obtained from scrapie-infected flocks from several loca- 
tions in Aragon (Spain); 7 of these exhibited clinical signs 
of scrapie while 4 were preclinical. Six animals were se- 
lected from flocks located in scrapie-free areas and were 
used as breed, genotype and age-matched controls. 

Sheep in the clinical disease stage were identified by 
assessing clinical signs associated with the disease [18]. 
Third eyelid [18] and rectal mucosa biopsies [19] were 
used to confirm this diagnosis and to identify animals in 
the preclinical disease stage. Postmortem examinations 
revealed no additional pathological findings. Mesenteric 
lymph node samples were divided into 2 halves; one was 
snap-frozen in liquid nitrogen and stored at -80°C until 
RNA extraction, and the other was formalin-fixed and 
paraffin-embedded for further histopathological analysis. 
Total RNA isolation, purification and quality control 
were performed as previously described [10,11]. 

Immunohistochemical detection of PrP Sc 

Immunohistochemical (IHC) analyses were performed 
using serial sections. Prion protein was detected, after for- 
mic acid treatment and proteinase K digestion, using the 
monoclonal primary antibody L42 (1:500; R-Biopharm), as 
previously described [20]. Negative controls were per- 
formed omitting the primary antibody from control and 
scrapie sections. 

The preparations were microscopically examined and 
global PrP Sc deposition was scored by quantifying the 
proportion of lymphoid follicles with PrP Sc deposits as 
follows: 1, <20%; 2, 20-50%; 3, >50% [19]. Significant dif- 
ferences between clinical, preclinical and control groups 
were identified using the Kruskal-Wallis test. 

Microarray hybridization 

The custom CVI 4x44K microarrays used in this study 
contained custom e Array- designed 60-mer probes of 
previously sequenced, normalized and subtracted cDNA 
libraries of ovine Peyer s Patch, obex and tonsil, supple- 
mented with the publicly available Ovis aries transcripts 
from the NCBI/EBI databases and the Agilent O. aries 
transcript catalog. All arrays were printed using Sureprint 
technology (Agilent Technologies). 

Preparation of the labeled cDNA probes and subsequent 
Genechip hybridizations were performed in accordance 
with the Agilent Technologies One- Color Microarray- 
Based Gene Expression Analysis guidelines, as previously 
described [10,11]. Hybridizations were performed in du- 
plicate, resulting in 14 microarrays for the clinical sam- 
ples, 8 microarrays for the preclinical samples and 12 for 
the negative control samples. Microarrays were scanned 
using the GenePix 4200AL Scanner (Axon Instruments) 
and GenePix Pro 6.0 software and the hybridization data 
were extracted using Agilent Feature Extraction software 
v9.5.3.1 (Agilent Technologies) before processing with 
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GeneSpring GX 10.0.2 (Agilent Technologies). Chip values 
were normalized using the 75th percentile value and ex- 
pression levels were calculated. Global mesenteric lymph 
node gene expression profiles for the clinical and preclin- 
ical disease stages were compared with the negative 
controls using a linear model that accounts for both tech- 
nical (random animal effects) and biological replicates 
(disease effects). A multiple testing correction developed 
by Benjamini and Hochberg was also applied. Only genes 
that exhibited at least a 2-fold change (FC) in expression 
between healthy and infected animals and had a P-value < 
0.05 were selected. These genes were clustered according 
to their Euclidean distance coefficient using PermutMatrix 
software [21]. A Megablast search of the GenBank data- 
base nr/nt was performed to identify genes that were simi- 
lar to the differentially expressed probes. The molecular 
functions of the genes were classified according to Gene 
Ontology (GO) using an on-line functional annotation 
tool (DAVID Bioinformatics Resources 6.7; NIAID/NIH, 
USA) [22,23]. 

Quantitative real-time PCR 

Quantitative real-time PCR (qRT-PCR) was performed to 
confirm the expression of 7 genes. Genes were selected 
based on previous studies describing their association with 
prion and other neurodegenerative diseases or their 
potential role in protein misfolding repair, the promotion 
of angiogenesis, or the systemic response of the organism 
to infection [24-30]. The PCR primer sequences used for 
the quantification of the genes encoding butyrobetaine 
(gamma), 2-oxoglutarate dioxygenase 1 (BBOX1), cerulo- 
plasmin (ferr oxidase) (CP), prefoldin subunit 2 (PFDN2), 
proteasome subunit, alpha type, 7 (PSMA7), serpin peptid- 
ase inhibitor, clade E (nexin, plasminogen activator inhibitor 
type 1), member 1 (SERPINE1), ubiquitin carboxyl-terminal 
esterase LI (ubiquitin thiolesterase) (UCHL1) and vascular 
endothelial growth factor A (VEGFA) are shown in Table 1. 
The RNA samples used for qRT-PCR were the same as 
those used for microarray experiments. The qRT-PCR 
assays were designed using Primer Express 2.0 software 
(Applied Biosystems) to select appropriate primer se- 
quences from known ovine or bovine sequences. When- 
ever possible, the exon-exon border was included to 
prevent amplification of genomic DNA in the PCR reac- 
tion. Complementary DNA (cDNA) was synthesized from 
1 \ig RNA using random hexamer primers with the Super- 
script First Standard Synthesis System for RT-PCR (Invi- 
trogen). To confirm the elimination of any remaining 
DNA, reverse transcription was performed with and with- 
out the enzyme. 

Quantitative RT-PCR was performed using SYBR® 
Green (PE Applied Biosystems) assays as previously de- 
scribed [10,11]. To improve the normalization accuracy, 
a normalization factor (NF) was used to determine the 
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Table 1 Genes analyzed by quantitative real-time PCR 



Gene 


Primer sequence 


Size 
(bp) 


Accession number 


BB0X1 


F: TGCAAACAATGTGGCTOCACA 
R: AAGCTGAACCCCAGGTGGAT 


85 


NM_ 


_001 101881.2** 


CP 


F: GCAGCCAGATACTGCAGGGAC 
R: CCCGCACTGGCTCACAGTATAT 


97 


NM_ 


_001 009733.1* 


PFDN2 


F: GCAGGTAATOCTGGOTCAA 
R: TTCAACTCCATCTCCAGCTCAG 


84 


NM_ 


_001080221.2** 


PSMA7 


F: TAATGTCTGCATGGCATOGC 
R: TGGCTCTGGCATOCACC 


81 


NM_ 


_001 034233.2** 


SERPINE1 


F: TGTACGTGTCGCAGGCGC 

R: TACAAGGGCTGTGGAGGAGGAC 


83 


NM_ 


_1 741 37.2** 


UCHL1 


F: AAOTGATGGACGGATGCOT 
R: TGCAGACOTGGCAGCGT 


84 


NM_ 


,0010461 72.2** 


VEGFA 


F: GGGCTGCTGTAATGACGAAAGT 
R: GGmGATCCGCATAATCTGCA 


81 


NM_ 


_001 025 110.1* 



*Ovine cDNA. 
**Bovine cDNA. 

Primers (F, forward; R, reverse) used for gene amplification, amplicon size and 
GenBank accession numbers for the bovine cDNA sequences used for primer 
design. 



expression level of each gene in each sample, as previ- 
ously described [31]. The NF was calculated using the 
geometric mean expression of 3 housekeeping genes 
(GAPDH, G6PDH and RPL32). The primers and PCR 
conditions used for the amplification of these house- 
keeping genes have been described previously [31,32]. 

The quantitative results obtained from the qRT-PCR 
assays were expressed as FC with respect to controls. 
Significant differences between groups were determined 
using a Students £-test (P < 0.05). 

The relationship between neuropathology and gene 
expression 

We used 2 different models to quantify the relationships 
between the PrP Sc deposition and gene expression. We 
first used a simple regression model, 

Ji = V + + et 

where y { is the gene expression data for the i th individual 
for each gene, \i is the general mean, b is slope of the re- 
gression analysis, p t is the PrP Sc deposition for the i th in- 
dividual and ei is the residual. In addition, we used 
another model, 

Jij ={*+Ti + b Pij + etj 

which includes a systematic effect associated with the 3 
categories (Preclinical, Clinical and Healthy). Here 
represents the gene expression for the j th individual at 
the i th category, Ti the effects of the i th treatment, p A j the 
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PrP Sc deposition and the residual The statistical sig- 
nificance of the slope associated with PrP Sc deposition 
was determined for both models. Analyses were per- 
formed using the R project for Statistical Computing. 

Results 

PrP Sc deposition in lymphoid tissue 

The accumulation of PrP Sc in the LRS was consistent 
with the features of classical scrapie. PrP Sc was detected 
within the primary and secondary lymphoid follicles of 
the lymph nodes. Some cells in the follicular mantle and 
paracortical area were also PrP Sc -positive. Interfollicular 
and sinusal cells were rarely labeled. As seen in other 
LRS tissues, distinct granular immunolabeling was de- 
tected in the cytoplasm of tingible body macrophages in 



both the dark and light zones of the lymphoid follicles. 
A granular deposition was observed in the vicinity of 
round cells, which were morphologically similar to im- 
mature B lymphocytes (Figure 1). None of the control 
animals displayed PrP Sc immunostaining and no statisti- 
cally significant differences were observed between the 2 
scrapie-infected groups (preclinical and clinical). 

Identification of differentially expressed genes in the 
mesenteric lymph node of sheep with natural scrapie 

A total of 234 probe sets displayed significant differences 
in expression of 2-fold or more between the control and 
scrapie (preclinical and clinical) groups. While the genes 
of Ovis aries are relatively poorly annotated, we used 
BLAST searches against publicly available databases to 



PrP Sc deposition 



C3 
> 
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Figure 1 Semi-quantitative assessment of PrP Sc deposition in mesenteric lymph nodes samples from control, preclinical and clinical 
scrapie sheep. Deposition was evaluated on a scale of 0 (negative) to 5 (maximum staining intensity). Significant differences were detected 
between control and scrapie groups using a Kruskal-Wallis test (**P<0.01) (A). PrP Sc deposition in the mesenteric lymph node from preclinical 
(B) and clinical (C) scrapie sheep. 
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identify a set of 177 known genes from the complete set 
of 234 differentially expressed probes. Some of these 234 
probes encoded the same gene. The microarray data 
were deposited in ArrayExpress and assigned the acces- 
sion number E-MTAB-1346. 

Based on the GO analysis of human homologs, 139 
genes had known functions (Table 2), of which 105 were 
differentially regulated in clinical scrapie-infected sheep 
(43 upregulated and 62 downregulated) and 50 in pre- 
clinical scrapie-infected sheep (12 upregulated and 38 
downregulated). The functional groups with the highest 
p values in the enrichment analysis are shown in Table 2; 
the most relevant terms were glycoprotein, extracellular 
region, disulfide bond and cell cycle. In addition, some of 
the annotated genes were grouped in 3 signaling path- 
ways: focal adhesion, PPAR signaling pathway and ECM- 
receptor interaction (Table 3). After hierarchical clustering 
using the expression levels of all annotated genes, we were 
able to discriminate the 3 experimental groups: clinical, 
preclinical and healthy (Figure 2). 

The relatively low number of animals analyzed and the 
fact that they were obtained from different geographical 
locations may limit the statistical power of our study. 
However, it is difficult to perform this kind of study with 
a larger number of animals due to the paucity of natural 
scrapie cases. Moreover, finding negative age-matched 
controls with the same PRNP genotype within the same 
scrapie- affected flocks is not possible since the disease 
attack-rate in these animals is close to 100%. 

Validation of gene expression profiling by quantitative 
RT-PCR 

To confirm the microarray results, the expression of 7 
genes (BBOX1, CP, PFDN2, PSMA7, SERPINE1, UCHL1 
and VEGFA) was analyzed by qRT-PCR. The expression 
of 4 of the selected targets (BBOX1, PFDN2, PSMA7 
and UCHL1) was altered in the same direction in both 
the clinical and preclinical stages of scrapie disease. Ex- 
pression levels of the other 3 genes (CP, SERPINE1 and 
VEGFA) were altered in opposite directions in the clin- 
ical and preclinical stages, respectively. In all cases, qRT- 
PCR confirmed the changes observed in the microarray 
hybridizations (Figure 3), revealing statistically signifi- 
cant differences between control and scrapie groups. 

Identification of prion deposition-related genes 

Two regression models were used to identify genes asso- 
ciated with prion deposition. One is a regression model 
that relates gene expression with prion deposition scores 
regardless of group (Preclinical, Clinical or Healthy), 
while the other includes a systematic effect associated 
with the 3 different animal categories and analyzes the 
remaining variation after subtracting the variability asso- 
ciated with each of the 3 groups. Natural conditions 



normally display high variability within groups. Using 
these regression analyses we identified genes whose ex- 
pression was associated with prion deposition, even if 
the changes in expression were not significant. We iden- 
tified 96 genes with known functions whose expression 
was significantly associated with PrP Sc deposition in 
both regression models (Figure 4). The slope value rep- 
resented was obtained using the second regression 
model. The gene ontology analysis revealed that genes 
associated with prion deposition were involved in DNA 
binding, the generation of precursor metabolites and en- 
ergy, ion binding, catabolic process and kinase functions. 
Only statistically significant slopes (p < 0.05) for both as- 
sociation models are presented in Figure 4. 

Discussion 

Mesenteric lymph nodes are one of the first organs in 
which PrP Sc accumulates [5]. To the best of our know- 
ledge, ours is the first study to use gene expression 
microarray technology to analyze mesenteric lymph 
nodes from naturally-infected sheep during the preclin- 
ical and clinical stages of scrapie. 

Our microarray hybridization analysis identified 139 
annotated genes of known functions that exhibited 
changes in expression of 2-fold or more with respect to 
controls during at least one of the 2 disease stages ana- 
lyzed. Of these 139 genes, 105 were differentially regu- 
lated in clinical scrapie-infected sheep and 50 in 
preclinical scrapie-infected sheep. These findings indi- 
cate that in the early disease stages, fewer genes are acti- 
vated or the changes in expression are not sufficiently 
large to be detected by microarray. PrP Sc accumulation 
in MLN peaked during the preclinical phase of scrapie 
(Figure 1), indicating that dysregulated genes that were 
not common to both disease stages could not be specif- 
ically linked to the prion accumulation process. 

Prion infection of neural cell lines has revealed a lim- 
ited transcriptional response and clear differences in 
transcriptional profiles between different cell lines and 
between cell cultures and whole brains from infected 
mice [35-37]. Similarly, our results differ from the pro- 
files described in the brains of sheep with natural scrapie 
[10,11], indicating the absence of universal changes asso- 
ciated with prion infection. 

The differentially expressed genes identified in the 
present study encode proteins located in almost all sub- 
cellular locations: membranes (glycoproteins), the extra- 
cellular matrix, the cytoplasm and chromosomes. These 
proteins participate in general cell processes such as the 
cell cycle, apoptosis, alternative splicing and chromo- 
some condensation. Genes associated with glycoprotein 
metabolism, the extracellular matrix and the cell cycle 
have been previously detected in leukocyte-depleted 
splenic cells from mice infected with scrapie [16]. We 



Table 2 GO terms for differentially expressed genes 



Term P value ProbeName FC (C vs NC) FC (PC vs NC) 



Glycoprotein 5.16E-04 


A_70_P042906 


2.22 


1.00 




CUST_5666_PI375351158 


-2.02 


-1.16 




A_70_P041077 


2.09 


1.39 




A_70_P030356 


1.48 


-2.06 




A_70_P051037 


-2.09 


2.85 




A_70_P036386 


-1.29 


-2.01 


Extracellular region 5.37E-04 


CUST_2502_PI375351 158 


2.26 


2.76 




A_70_P006426 


4.26 


1.06 




A_70_P069571 


2.05 


1.93 




A_70_P038041 


1.95 


-3.63 




CUST_3439_PI375351158 


2.70 


1.38 




A_70_P035871 


1.03 


-3.29 




A_70_P039451 


-1.57 


-2.28 




A_70_P050922 


2.19 


-1.06 


Disulfide bond 1.10E-03 


CUST_598_PI375351 158 


2.28 


1.33 




A_70_P070736 


-3.83 


-1.06 


Cell cycle 1.17E-03 


A_70_P058586 


-2.85 


-1.55 




A_70_P057841 


-1.61 


-2.51 




A_70_P053686 


1.11 


-2.01 




A_70_P060136 


-1.48 


-2.04 




A_70_P027076 


-2.79 


-1.28 




A_70_P007666 


-2.87 


-1.61 




A_70_P01 1891 


-2.34 


-1.29 




A_/U_rUUyjJ 1 


1 1 3 

— Z. I 0 


1 H7 

— I .u/ 




A_70_P037556 


-2.12 


1.01 




A_70_P034256 


-2.20 


1.05 




A_70_P061801 


-2.01 


-1.39 


Condensed chromosome 1.71E-03 


CUST_7895_PI375351 158 


-2.79 


-1.15 




A_70_P003816 


-2.20 


-1.67 



Gene symbol Gene name Citation 



AQP1 Aquaporin 1 (Colton blood group) 

AQP4 Aquaporin 4 [10] 

PTGER3 Prostaglandin E receptor 3 (subtype EP3) 

SLC1 1 A1 Solute carrier family 1 1 (proton-coupled divalent 
metal ion transporters), member 1 

SLC2A5 Solute carrier family 2 (facilitated glucose/fructose [10] 
transporter), member 5 

UCHL1* Ubiquitin carboxyl-terminal esterase L1 (ubiquitin thiolesterase) 

BMPER BMP binding endothelial regulator 

FGL1 Fibrinogen-like 1 

FIBIN Fin bud initiation factor homolog (zebrafish) 

HP Haptoglobin [10] 

PIGR Polymeric immunoglobulin receptor 

SERPINE1* Serpin peptidase inhibitor, clade E (nexin, plasminogen 
activator inhibitor type 1), member 1 

SFN Stratifin 

VEGFA* Vascular endothelial growth factor A 

FABP5 Fatty acid binding protein 5 (psoriasis-associated) 

TXNRD1 Thioredoxin reductase 1 

PSMA7* Proteasome (prosome, macropain) subunit, alpha type, 7 

BRCA2 Breast cancer 2, early onset 

CCNDBP1 Cyclin D-type binding-protein 1 

CGREF1 Cell growth regulator with EF-hand domain 1 

EID1 EP300 interacting inhibitor of differentiation 1 

KIF2C Kinesin family member 2C 

MND1 Meiotic nuclear divisions 1 homolog (S. cerevisiae) 

NCAPG Non-SMC condensin I complex, subunit G 

NDC80 NDC80 homolog, kinetochore complex 
component (S. cerevisiae) 

TFDP2 Transcription factor Dp-2 (E2F dimerization partner 2) 

UBE2C Ubiquitin-conjugating enzyme E2C 

HMGB1 High-mobility group box 1 

TEX12 Testis expressed 1 2 



Table 2 GO terms for differentially expressed genes (Continued) 



Alternative splicing 


2.04E-03 


A_70_P059171 


-2.57 


-2.17 






A_70_P033441 


1.13 


2.27 






A_70_P039536 


-2.86 


-1.34 






A_70_P060791 


2.25 


1.38 






A_70_P015446 


-1.52 


-2.13 


Chromosome 


2.20E-03 


A_70_P065046 


-2.68 


-1.32 






A_70_P040916 


-2.62 


-1.67 






CUST_12205_PI375351158 


-2.04 


-1.20 


Extracellular matrix 


3.15E-03 


A_70_P021711 


-2.19 


-1.69 






CUST_1 1 292_PI37535 1 1 58 


2.40 


1.12 






A_70_P007356 


2.20 


1.26 






CUST_7449_PI375351 158 


1.46 


-2.05 






A_70_P032476 


2.45 


-1.35 






A_70_P058626 


2.27 


1.13 


Signal peptide 


3.25E-03 


A_70_P064081 


2.44 


-2.18 






A_70_P061221 


1.20 


-2.43 






A_70_P040831 


-2.09 


-1.07 


Tumor suppressor 


1 .47E-02 


A_70_P035971 


2.05 


-1.75 






A_70_P026737 


3.11 


-1.41 






A_70_P061691 


1.26 


2.26 


Cytoplasm 


1 .82E-02 


A_70_P016337 


-1.78 


-2.72 






CUST_12176_PI375351 158 


-1.04 


-2.92 






A_70_P019656 


1.01 


3.59 






A_70_P045301 


3.70 


2.33 






A_70_P069971 


-2.09 


-1.73 






CUST_5823_PI375351 158 


1.12 


2.07 






A_70_P0 18021 


-3.25 


-1.15 






A_70_P054126 


-2.16 


-1.11 






A_70_P044246 


1.26 


-2.14 






A_70_P055106 


2.05 


1.07 






A_70_P025726 


-2.39 


-1.60 






A_70_P060991 


-2.80 


1.21 



CLTA Clathrin, light chain (Lea) 

MEST Mesoderm specific transcript homolog (mouse) 

MX1 Myxovirus (influenza virus) resistance 1, 
interferon-inducible protein p78 (mouse) 

PEG3 Paternally expressed 3 

UNC50 Unc-50 homolog (C. elegans) 

JRKL Jerky homolog-like (mouse) 

ORC1L Origin recognition complex, subunit 1-like (yeast) 

TOP2A Topoisomerase (DNA) II alpha 170 kDa 

COL1A1 Collagen, type I, alpha 1 

CRISPLD2 Cysteine-rich secretory protein LCCL domain containing 2 

DCN Decorin 

DPT Dermatopontin 

MFAP5 Microfibrillar associated protein 5 

VWF Von Willebrand factor 

CA4 Carbonic anhydrase IV 

P4HA1 Prolyl 4-hydroxylase, alpha polypeptide I 

SELL Selectin L 

VWA5A Von Willebrand factor A domain containing 5A 

WT1 Wilms tumor 1 

HPGD Hydroxyprostaglandin dehydrogenase 15-(NAD) 

HMMR Hyaluronan-mediated motility receptor (RHAMM) 

ACTG2 Actin, gamma 2, smooth muscle, enteric 

ADH1A Alcohol dehydrogenase 1A (class I), alpha polypeptide 

BBOX1* Butyrobetaine (gamma), 2-oxoglutarate dioxygenase 
(gamma-butyrobetaine hydroxylase) 1 

CDC42EP5 CDC42 effector protein (Rho GTPase binding) 5 

NIN Ninein (GSK3B interacting protein) 

PAIP2 Poly(A) binding protein interacting protein 2 

PHAX Phosphorylated adaptor for RNA export 

RRAGD Ras-related GTP binding D 

SORBS2 Sorbin and SH3 domain containing 2 

STMN2 Stathmin-like 2 

CA3 Carbonic anhydrase III, muscle specific 



[33] 
[10,16] 

[10,16] 

[10] 

[10,11] 
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Table 2 GO terms for differentially expressed genes (Continued) 



Fatty acid metabolism 


2.26E-02 


CUST_8730_PI375351 158 


1.06 


-2.05 






A_70_P039986 


-1.27 


2.39 


Phosphatidic acid 


2.86E-02 


A_70_P039726 


-2.50 


-1.10 


metabolic process 














A_70_P025121 


-3.34 


-1.30 


Protein dimerization activity 


2.91 E-02 


CUST_206_PI396851 701 


2.04 


-2.28 






CUST_10535_PI375351158 


-2.41 


-2.41 


Zinc finger, C2H2-type 


3.09E-02 


CUST_13307_PI375351158 


-2.34 


-1.26 






A_70_P036056 


-2.65 


-1.71 






A_70_P065371 


-2.27 


-1.13 






A_70_P0 14286 


1.33 


2.01 


Cytoplasmic 


3.21 E-02 


A_70_P062831 


-3.71 


-1.48 


membrane-bounded 




A_70_P044621 


2.27 


-1.04 


vesicle 




Metal ion transport 


3.39E-02 


A_70_P030126 


-2.36 


1.24 






CUST_12632_PI375351158 


2.28 


1.50 






A_70_P027206 


2.15 


2.10 






A_70_P050456 


-2.65 


-1.58 


Antimicrobial 


5.05E-02 


CUST_139_PI396851701 


15.26 


-1.74 






CUST_3993_PI375351158 


5.96 


-1.25 


Apoptosis 


5.45E-02 


CUST_260_PI396851 701 


-2.04 


-1.55 






A_70_P022526 


-2.29 


-1.16 






A_70_P049526 


-1.12 


-2.70 






A_70_P069836 


-3.14 


-1.31 






A_70_P067091 


-2.90 


-1.55 


Immunoglobulin-like fold 


5.55E-02 


A_70_P036361 


2.06 


2.09 






CUST_9656_PI375351 158 


-2.03 


1.11 






CUST_12141_PI375351158 


2.64 


-1.00 






A_70_P0 16986 


-2.23 


1.16 


Secreted 


5.70E-02 


A 70 P04R761 


1 .34 


-2.27 






CUST_9580_PI375351158 


1.13 


-2.49 






CUST_10101_PI375351158 


-1.07 


-2.53 






A_70_P034356 


2.88 


1.19 






A_70_P063891 


2.38 


-1.22 



Acyl-CoA synthetase long-chain family member 4 

Solute carrier family 27 (fatty acid transporter), member 2 

1-acylglycerol-3-phosphate O-acyltransferase 1 
(lysophosphatidic acid acyltransferase, alpha) 

Phospholipase A2, group IIA (platelets, synovial fluid) 

v-fos FBJ murine osteosarcoma viral oncogene homolog 

Nuclear factor (erythroid-derived 2)-like 1 

ATM interactor 

Kruppel-like factor 6 

Zinc finger, CCHC domain containing 6 

Zinc finger protein 45 

EH-domain containing 3 

Mai, T-cell differentiation protein-like 

Ceruloplasmin (ferroxidase) 

Potassium channel tetramerisation domain containing 12 
Solute carrier family 30 (zinc transporter), member 1 
Solute carrier family 38, member 10 
Cathelicidin 1 

S100 calcium binding protein A8 

Caspase recruitment domain family, member 6 

Optic atrophy 1 (autosomal dominant) 

PERP, TP53 apoptosis effector 

Ring finger protein 130 

Ubiquitin-conjugating enzyme E2Z 

CD1e molecule 

Fibroblast growth factor receptor 2 

Ig kappa chain 

T cell receptor delta locus 

Serine peptidase inhibitor, Kazal type 5 

Multimerin 1 

Nephronectin 

Sema domain, immunoglobulin domain (Ig), short basic 
domain, secreted, (semaphorin) 3G 

Procollagen C-endopeptidase enhancer 2 



Table 2 GO terms for differentially expressed genes (Continued) 



GTP binding 7.49E-02 A_70_P019676 4.33 1.16 

A_70_P046951 -2.24 -1.02 

A_70_P008796 -2.06 1.15 

A_70_P000126 -3.63 -1.82 

Disulfide bond 7.58E-02 A_70_P018376 2.56 -1.08 

A_70_P041321 -2.15 -2.57 

Structural molecule activity 8.04E-02 A_70_P028306 -2.03 -1.05 

CUST_13373_PI375351158 -2.15 -1.29 

CUST_7057_PI375351158 -2.33 -1.07 

A_70_P031441 -2.43 -1.79 

Positive regulation of 9.48E-02 A_70_P058551 8.25 -1.12 

inflammatory response CUST_9779_PI375351 158 -1.40 -2.38 

Other CUST_1 3467_PI37535 1 1 58 1.29 -2.14 

A_70_P055306 2.77 1.50 

A_70_P049406 -1.99 -2.32 

A_70_P070931 -2.87 -1.18 

A_70_P036101 1.22 -3.18 

A_70_P01 1091 -2.74 -2.01 

A_70_P039656 17.20 -1.52 

A_70_P061341 -1.03 -2.94 

A_70_P061266 2.36 1.03 

A_70_P051377 14.99 -1.48 

A_70_P015101 -2.70 -1.03 

A_70_P012541 -1.02 -2.36 

A_70_P002841 2.03 1.48 

A_70_P008841 -2.18 -1.23 

A_70_P060756 -2.12 1.08 

A_70_P062951 2.05 1.12 

A_70_P051181 2.76 1.28 

CUST_46_PI396851701 -2.06 -2.15 

A_70_P024671 2.46 1.31 



Acyl-CoA synthetase medium-chain family member 1 

ADP-ribosylation factor-like 4C 

Dedicator of cytokinesis 1 1 

Guanylate binding protein 4 

Fibroblast activation protein, alpha 

Immunoglobulin heavy constant alpha 2 (A2m marker) 

Lamin B1 

Ribosomal protein 17-like; similar to 60S ribosomal 
Protein L17 (L23); similar to ribosomal protein LI 7; 
similar to Rpl17 protein; similar to mCG3798 

Similar to Ribosomal_L22 domain containing protein RGD1 359290 

Vinculin 

Fatty acid binding protein 4, adipocyte 

Indoleamine 2,3-dioxygenase 1 

Fructose-1,6-bisphosphatase 2 

Solute carrier family 43, member 1 

Eukaryotic translation initiation factor 5 

Nudix (nucleoside diphosphate linked moiety X)-type motif 19 

NOP10 ribonucleoprotein homolog (yeast) 

Prefoldin subunit 2 

5 kDa bactinecin precursor 

Chromosome 11 open reading frame 10 

Chromosome 5 open reading frame 43 

Procyclic dodecapeptide 

Coiled-coil domain containing 82 

Calponin 1, basic, smooth muscle 

Endothelial PAS domain protein 1 

Epithelial stromal interaction 1 (breast) 

KIAA1128 

Fibronectin type III domain containing 3B 
NAD-dependent 1 1-beta-hydroxysteroid dehydrogenase 
Immunoglobulin heavy constant alpha 
Inhibin, beta B 



Table 2 GO terms for differentially expressed genes (Continued) 



A_70_P033656 


-2.34 


-1.77 


ISG 12(A) 


ISG1 2(a) protein-like 




A_70_P012537 


-2.19 


-3.58 


LGALS15 


Galectin 15 




CUST_3620_PI375351158 


-1.58 


-2.22 


LGALS16 


Beta-galactoside-binding lectin 




A_70_P020131 


-1.10 


-3.84 


LGALS9C 


Lectin, galactoside-binding, soluble, 9C 




CUST_6089_PI375351 158 


-2.02 


-1.71 


LOC100140018 


Similar to cytochrome P450, family 2, subfamily J 




A_70_P016801 


-2.20 


-1.48 


LOC443321 


Lysozyme 2a precursor 




A_70_P026671 


-1.91 


-2.05 


LOC524810 


IgM 




A_70_P023026 


3.09 


1.74 


LOC6 15697 


Similar to cytochrome P450 




A_70_P039266 


-2.03 


-1.83 


LOC654331 


Pancreatitis-associated protein I 


[10] 


A_70_P039271 


-3.94 


-7.40 


MCP-4 


Mast cell proteinase-4 




CUST_1 1 558_PI37535 1 1 58 


2.00 


2.23 


MGC 140754 


Hypothetical LOC508613 




CUST_9725_PI375351158 


-2.26 


-3.77 


PLAC8 


Placenta-specific 8 




CUST_10114_PI375351158 


2.56 


1.15 


RN18S1 


18S ribosomal RNA 




A_70_P030132 


-2.40 


-1.71 


RTP4 


Receptor (chemosensory) transporter protein 4 


[34] 


A_70_P030426 


10.95 


1.39 


SC5 


Cathelin-related prepropeptide 





*Genes selected for validation by quantitative RT-PCR. 

BLAST results for clones exhibiting significant alterations in gene expression (P<0.05 and > 2-fold change) in preclinical and clinical stages of natural scrapie. The GO was determined using DAVID Bioinformatics 
Resources 6.7 (NIAID/NIH, USA), an on-line functional annotation tool. Only genes with a known GO are shown. (C, Clinical; FC, Fold change; NC, Negative control; PC, Preclinical). Previous studies describing altered 
expression of the differentially expressed genes reported here in the brain and lymphoid tissues of experimentally-infected sheep are cited. 
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Table 3 Signaling pathways associated with differentially expressed genes 



Term 


P value 


ProbeName 


FC (C vs NC) 


FC (PC vs NC) 


Gene 


Gene name 
















symbol 




DDAD r i r-\ ~i 1 i r> 

rrAK signaling 


I .oot-Uz 


CUST 


_598_PI375351 158 


Z.ZO 


1 o o 

1 .55 


r ADnr 
rAbrj 


Fatty acid binding protein 5 (psoriasis-associated) 


pathway 




CUST_8730_PI375351158 


I .(Jo 


— z.Uj 


ALoL4 


Acyl-CoA synthetase long-chain family member 4 






A_ 


_70_ 


_P058551 


Q 1C 
O.ZJ 


— I . I z 


C ADD/1 

rAbr4 


Fatty acid binding protein 4, adipocyte 






A_ 


_70_ 


_P039986 


— I .Z/ 


z.oy 


oLLz/Az 


Solute carrier family 27 (fatty acid transporter), 
member 2 


ECM-receptor 


2.53E-02 


A_ 


_70_ 


_P021711 


-2.19 


-1.69 


C0L1A1 


Collagen, type I, alpha 1 


interaction 




A_ 


_70_ 


_P058626 


2.27 


1.13 


VWF 


von Willebrand factor 






A_ 


_70_ 


_P016337 


-1.78 


-2.72 


HMMR 


Hyaluronan-mediated motility receptor (RHAMM) 


Focal adhesion 


6.73E-02 


A_ 


_70_ 


_P021711 


-2.19 


-1.69 


C0L1A1 


Collagen, type I, alpha 1 






A_ 


_70_ 


.P058626 


2.27 


1.13 


VWF 


von Willebrand factor 






A_ 


_70_ 


_P050922 


2.19 


-1.06 


VEGFA 


Vascular endothelial growth factor A 






A_ 


_70_ 


_P031441 


-2.43 


-1.79 


VCL 


Vinculin 



Signaling pathways with which differentially expressed genes (P<0.05 and > 2-fold change) are associated were identified by enrichment analysis using the 
on-line functional annotation (DAVID Bioinformatics Resources 6.7; NIAID/NIH, USA). (C, Clinical; FC, Fold change; NC, Negative control; PC, Preclinical). 



found that gene encoding decorin (DCN), a protein lo- 
cated in the extracellular matrix, was overexpressed in 
scrapie animals, in line with previous reports of DCN 
upregulation at both the gene and protein level in cells 
isolated from scrapie-infected mouse spleen [16]. Our 
results confirm the differential regulation of DCN in the 
MLN of sheep with natural scrapie. Some of differen- 
tially transcribed genes that we identified encode pro- 
teins that may be involved in the pathogenesis of protein 
misfolding diseases. Examples include proteins involved 
in disulfide bonding, GTP binding and protein dimeriza- 
tion. Differential expression of some of these genes has 
been reported in other tissues in scrapie-infected mice 
[12,14,38] and in natural scrapie [10,11]. 

The GO enrichment analysis of our microarray data 
identified 3 signaling pathways that may be involved in 
the pathogenesis of scrapie; extracellular matrix (ECM)- 
receptor interaction, focal adhesion and the peroxisome 
proliferator- activated receptor (PPAR) signaling pathway. 
We and others have previously reported altered expres- 
sion of COL1A2, COL3A1 and COL12A1 and other 
ECM related genes in the CNS of sheep infected with 
scrapie and other prion diseases [10,11,39,40] and in 
leukocyte-depleted splenic cells from scrapie-infected 
mice [16]. Moreover, prion protein regulates pi integrin 
signaling activity in prion-infected neurons [41] and the 
expression of genes of the PPAR signaling pathway is al- 
tered in Alzheimer's disease (AD) [42,43]. Our findings 
indicate alterations in the extracellular matrix, focal ad- 
hesion and in fatty acid metabolism via the PPAR signal- 
ing pathway in non-neural tissues susceptible to prion 
infection. 

In addition to analyzing differential gene expression, 
we investigated the effect of prion deposition on LRS 
gene expression in order to further our understanding of 



the prion-specific pathogenic pathways in scrapie. Two 
regression models (see Methods section) were used to 
identify possible associations between gene expression 
data and the PrP Sc deposition score, as determined by 
immunohistochemistry. These statistical association ap- 
proaches can identify genes that may be implicated in 
prion-specific processes, although the differences in their 
expression levels between control and infected animals 
were not significant. Except for CUT A, none of the genes 
identified in the association analysis exhibited significant 
differences in expression in the array data. The gene 
ontology analysis revealed that genes associated with 
prion deposition are involved in DNA binding, the gen- 
eration of precursor metabolites and energy, ion binding, 
catabolic processes and kinase functions. Among the 
genes that were closely associated with prion protein de- 
position were cathepsin K (CTSK) and clathrin, light 
chain A (CLTA). Several members of the cathepsin family 
(including cathepsins C, D, H, S and Z) are differentially 
expressed in prion-affected brain tissues [12,14,38,44,45]. 
We observed a negative association between CTSK expres- 
sion and prion deposition, suggesting a possible proteolytic 
effect of cathepsin K on prion protein. Furthermore, CLTA 
expression was negatively associated with prion depos- 
ition. The expression of this gene is also altered in many 
neurodegenerative diseases such as Alzheimer's and 
Huntington's diseases [46,47], and clathrin-mediated path- 
ways are intimately involved in amyloid formation in 
neurodegenerative diseases [48]. These observations sug- 
gest that CLTA may also be involved in prion amyloid 
formation. 

Seven of the differentially expressed genes detected 
were selected for further analysis by qRT-PCR. These 
genes were chosen based on their involvement in protein 
misfolding pathways (PFDN2, UCHL1 and PSMA7), the 
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Figure 2 Condition trees generated by cluster analysis. A 

hierarchical cluster analysis (Euclidean distance clustering algorithm) 

was performed using PermutMatrix [21], which identified 139 genes 

whose expression differed significantly (> 2 fold change with 

respect to controls) in at least at one of the 2 disease stages 

analyzed (preclinical and clinical). Colors represent the level of 

expression. Sample information is listed across the top. The names 

of the known genes are indicated. Note the distinct patterns of 

gene expression in control, preclinical and clinical animals. 
I J 

regulation of angiogenesis (VEGFA and SERPINE1), iron 
metabolism and neurotoxicity (CP) and antioxidant 
mechanisms (BB0X1). Comparable fold changes in the 
expression of all 7 genes, in both preclinical and clinical 
disease stages, were observed in the microarray and 
qRT-PCR analyses, confirming the validity of both the 
array data and the alignment analysis. 

Here, we will briefly discuss how these genes may be 
involved in PrP Sc formation, dissemination and toxicity 
and in the response of the organism to prion infection. 
Genes encoding proteins involved in molecular pathways 
that regulate protein misfolding (PSMA7, UCHL1 and 
PFDN2) were downregulated in scrapie-infected animals. 
Furthermore, our association study revealed a positive 
association between proteasome subunit, beta type, 8 
(PSMB8) and PrP Sc deposition. The proteins encoded by 
PSMA7, UCHL1 and PSMB8 form part of the ubiquitin 
proteasome system (UPS), a complex that regulates the 
degradation of incomplete, damaged or misfolded pro- 
teins [49,50]. The UCHL1 protein has been implicated 
in protein misfolding diseases. UCHL1 is abundantly 
expressed in nervous system tissues [51], and gene vari- 
ants and changes in the activity of UCHL1 have been as- 
sociated with neurodegenerative disorders such as 
Alzheimer's, Parkinsons and Huntington's diseases 
[27,52-55]. In prion diseases, the functional capacity of 
the UPS is impaired by the direct interaction of (3-sheet- 
rich PrP with the 20S core particle of the proteasome 
complex, inhibiting substrate entry [56]. UPS impair- 
ment may enable the conversion of cytosolic PrP c to an 
abnormal, "PrP Sc -like" form [57]. Increased levels of ubi- 
quitin conjugates are found in scrapie-infected mouse 
brains, and are correlated with decreased proteasome 
function [58]. The downregulation of PSMA7 and 
UCHL1 genes and the positive association of PrP Sc de- 
position with PSMB8 expression in both preclinical and 
clinical stages of natural scrapie in the LRS of sheep is in 
line with the impaired proteasome function described in 
prion diseases of the nervous system [28]. Moreover, 
prefoldin (PFD) is likely to bind to substrate proteins 
that exist in an unfolded state and to transfer these pro- 
teins to the cytosolic chaperonin-containing TCP-1 
(CCT) for functional folding [59]. The role of this pro- 
tein in protein misfolding diseases has been investigated 
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Fold Change during preclinical stage Fold Change during clinical stage 

Figure 3 Confirmation of microarray results by quantitative real-time PCR. Differential expression of selected genes analyzed by microarray 
and quantitative RT-PCR in mesenteric lymph node samples from preclinical and clinical scrapie sheep: gamma-butyrobetaine hydroxylase 
(BBOX1), ceruloplasmin (CP), prefoldin subunit 2 (PFDN2), proteasome subunit alpha type-7 (PSMA7), plasminogen activator inhibitor-1 (SERPINE1), 
ubiquitin carboxy-terminal hydrolase L1 (UCHL1), and vascular endothelial growth factor (VEGFA). 

K J 



and prefoldin subunit 2 (PFDN2) has been implicated in 
AD [60]. Our findings reveal downregulation of PFDN2 
expression in both preclinical and clinical stages of nat- 
ural scrapie in the LRS of sheep. The observed downreg- 
ulation of PSMA7 and UCHL1 may facilitate prion 
conversion, implicating these genes in the pathogenic 
mechanism of the scrapie. 

The expression of SERPINE 1 or plasminogen activator 
inhibitor-1 (PAI-1) was downregulated in preclinical scra- 
pie animals and upregulated in the clinical phase of the 
disease, suggesting a dual role for this gene in scrapie. 
SERPINE 1 is an inhibitor of 2 types of plasminogen acti- 
vators; tissue-type plasminogen (tPA) and urokinase-type 
plasminogen (uPA) activators. Plasminogen was recently 
show to stimulate PrP misfolding by accelerating the 
conversion of PrP c to PrP Sc [26]. The downregulation of 
SERPINE 1 in preclinical scrapie may inhibit PrP Sc forma- 
tion in the early stages of the disease. SERPINE 1 is also in- 
volved in angiogenesis, as evidenced by the defective 
angiogenesis seen in tumor-bearing SERPINE 1 -deficient 
mice [61]. Accordingly, decreased expression of SERPINE1 
may represent a protective response to PrP Sc propagation 
during the preclinical stage of natural scrapie. However, 
we observed a slight increase in SERPINE1 expression in 
clinical scrapie-infected sheep as compared with controls 
(FC = 1.25). The expression of vascular endothelial growth 
factor (VEGF, also known as VEGFA), which is also in- 
volved in angiogenesis, was altered in scrapie. This gene 
mediates a critical rate-limiting step in physiological 
angiogenesis and plays a key role in tumor progression 
and angiogenesis and its expression by primary tumor 
promotes lymphangiogenesis and facilitates metastasis to 
sentinel lymph nodes [62]. As observed for SERPINE1, we 



detected a slight decrease in VEGFA expression in lymph- 
atic tissues in the preclinical stage of natural scrapie, 
followed by a moderate increase during the clinical stage 
of the disease. Lymphatic vessels are one of the proposed 
routes of prion dissemination [63]. Accordingly, the differ- 
ential regulation of SERPINE1 and VEGFA observed in 
scrapie may reflect a role in prion propagation. However, 
further studies will be necessary to clarify the role of these 
proteins in prion dissemination. 

In addition to genes involved in prion misfolding and 
propagation, we analyzed the expression of the cerulo- 
plasmin (CP) gene, which is involved in neurotoxicity. 
CP is a ferroxidase that contains 95% of the copper 
present in blood plasma, and is thus a key mediator of 
copper transport and metabolism [64]. Copper plays an 
important role in prion integrity; in vitro, copper en- 
hances renaturation and stabilization of PrP Sc , restoring 
protease resistance and infectivity [64]. Increased CP ex- 
pression has been reported in leukocyte- depleted splenic 
cells from clinical scrapie-infected mice [16]. However, 
we found that CP expression was significantly downreg- 
ulated in the lymph nodes of scrapie-infected animals in 
the clinical stage of the disease. In patients with acerulo- 
plasminemia, the predominant clinical symptoms are 
neurologic and include dystonia, dysarthria, and subcor- 
tical dementia due to progressive degeneration of the 
basal ganglia, as well as selective neuronal loss and in- 
creased iron content in the microglia, neurons and basal 
ganglia [29]. Ceruloplasmin also plays an essential role 
in normal brain iron metabolism and protects against 
neuronal loss [65]. Our findings suggest that CP defi- 
ciency may contribute to the neuronal degeneration ob- 
served in scrapie, although further analyses of brain 
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P1 


P2 


Slope of regression I 


Actin cytoskeleton 


CTTNBP2NL 


CTTNBP2 N-teri 


1.50E-02 


3,06E-03 


^27^ 




MYL12B 


Myosin, light chain 12B, regulatory 


1.30E-02 


2.09E-02 


-2,26 ■ 




TBCC 


Tubulin folding cofactor C 


1 ,48E-02 


3.18E-02 


■ 2,03 


Cellular macromolecular complex 


H2AFY 


Core histone macro-H2A 


1 ,65E-03 


1 ,52E-02 


-2,43 


assembly 


TSPYL4 


TSPY-like 4 


1 ,25E-02 


3.24E-02 


-2,02 ■ 


Cellular macromolecule catabolic process 


ANAPC13 


Anaphase promoting complex subunit 13 


2.87E-02 


4,89E-02 


-1,78 ■ 




ERCC2 


Excision repair cross-complementing, group 2 


4,78E-02 


8,36E-03 


2,74 




KCMF1 


Potassium channel modulatory factor 1 


1 ,37E-02 


1 ,88E-02 


-2,32 ^| 




PSMB8 


Proteasome subunit, beta type, 8 


4,09E-03 


2.16E-02 


■ 2,24 




WWP2 


WW domain containing E3 ubiquitin protein ligase 2 


2,27E-03 


1,81 E-02 


■ 2,34 


Cytoskeletal part 


LATS1 


LATS, large tumor suppressor, homolog 1 (Drosophila) 


3,24E-02 


1.79E-02 


■ 2,34 




NME2 


NME/NM23 nucleoside diphosphate kinase 2 


1.66E-03 


3,32E-03 


-3,23 




RPL22 


Ribosomal protein L22 


4,81 E-02 


3,93E-02 


-1,91 ■ 


DNA binding 


ATF6B 


Activating transcription factor 6 beta 


4.95E-02 


7,44E-03 


-2,81 ■ 




NFYA 


Nuclear transcription factor Y, alpha 


3,44E-02 


4,80E-02 


-1,79 ■ 




PDZD11 


PDZ domain containing 11 


6.97E-03 


1.12E-02 


-2,59 ■ 




PUF60 


Poly-U binding splicing factor 60KDa 


3,98E-02 


3.73E-02 


-1,94 ■ 




R0RA 


RAR-related orphan receptor A 


4,74E-02 


3,84E-02 


-1,92 ■ 




TSC22D2 


TSC22 domain family, member 2 


4.56E-02 


3.28E-02 


-2,01 ■ 


Extracellular region 


B2M 


Ovis aries beta-2-microglobulin 


2.28E-02 


3,85E-02 


■ 1,92 


Extracellular region 


PDE3B 


Phosphodiesterase 3B, cGMP-inhibited 


1.78E-02 


1 ,79E-02 


-2,34 ^| 


Generation of precursor metabolites and 


ATP5G1 


ATP synthase, subunit C1 (subunit 9) 


1 ,28E-02 


1 ,90E-02 


-2,31 ■ 


energy 


NDUFB5 


NADH dehydrogenase (ubiquinone) 1 beta subcomplex, 5 


2,00E-02 


1 ,34E-02 


-2,49 ^| 




PDHA1 


Pyruvate dehydrogenase (lipoamide) alpha 1 


2,64E-02 


2,94E-02 


■ 2,07 




TPI1 


Triosephosphate isomerase 1 


1 ,38E-02 


1 ,77E-02 


2 ,35 


Integral to membrane 


FAM73B 


Family with sequence similarity 73, member B 


2,31 E-02 


4,51 E-02 


-1,83 ■ 


Ion binding 


HARBI1 


Harbinger transposase derived 1 


2,57E-02 


1 ,35E-02 


-2,49 ■ 




SNRPC 


Small nuclear ribonucleoprotein polypeptide C 


7.02E-03 


5,39E-03 


-2,97 




SRI 


Sorcin 


1.55E-02 


4,61 E-03 


3,05 




ZBTB34 


Zinc finger and BTB domain containing 34 


2,34E-02 


3,81 E-02 


■ 1,93 




ZC3HC1 


Zinc finger, C3HC-type containing 1 


8,47E-03 


4,41 E-02 


■ 1,84 


Kinase 


EEF2K 


Eukaryotic elongation factor-2 kinase 


1 ,50E-02 


4,53E-02 


-1,83 ■ 




PFKP 


Phosphofructokinase, platelet 


2,71 E-02 


2,61 E-02 


■ 2,14 




PHKA1 


Phosphorylase kinase, alpha 1 (muscle) 


3,13E-02 


3,37E-02 


■ 2,00 




SLK 


STE20-like kinase (yeast) 


3,63E-02 


4,25E-02 


-1,87 ■ 


Lymphocyte activation 


CDK2AP2 


Cyclin-dependent kinase 2 associated protein 2 


2,78E-02 


2.16E-02 


■ 2,24 




CXCL12 


Chemokine (C-X-C motif) ligand 12 


1,91 E-02 


1 ,99E-02 


^| 2,28 


Metal ion binding 


ASXL2 


Similar to additional sex combs like 2 


2,93E-02 


4,92E-03 


-3,02 




FEM1B 


Fem-1 homolog b (C. elegans) 


2,85E-03 


3,25E-03 


3,24 


Metal-binding 


GRAMD1C 


GRAM domain containing 1C 


9,71 E-04 


2.52E-02 


■ 2,16 


Mitochondrion 


CHCHD3 


Coiled-coil-helix-coiled-coil-helix domain containing 3 


4,25E-03 


1 ,95E-02 


2,30 




DUT 


Deoxyuridine triphosphatase 


2,20E-02 


1.14E-02 


-2,58 




MRPL22 


Mitochondrial ribosomal protein L22 


3,15E-02 


8,44E-03 


2,74 




TIMM10 


Translocase of inner mitochondrial membrane 10 


1 ,54E-02 


2,46E-02 


-2,17 ^| 


Non-membrane-bounded organelle 


CENPT 


Centromere protein T 


1 ,84E-02 


2,65E-02 


-2,13 ■ 




GAPT 


GRB2-binding adaptor protein, transmembrane 


1.07E-02 


3,88E-02 


-1,92 ■ 




ZWINT 


ZW10 interactor 


2,97E-02 


3,24E-02 


-2,02 B 


Nucleotide binding 


GMPS 


Guanine monphosphate synthetase 


4,90E-02 


4.87E-02 


-1,79 ■ 




GSPT1 


G1 to S phase transition 1 


2,62E-02 


3,21 E-02 


■ 2,02 




RECQL4 


RecQ protein-like 4 


2.97E-02 


1 ,04E-02 


-2,63 


Organelle envelope 


CLTA 


Clathrin, light chain A 


1 ,42E-03 


3,84E-02 


-1,92 ■ 




XP07 


Exportin 7 


2,51 E-02 


1.32E-02 


2,51 


Phosphorus metabolic process 


DUSP10 


Dual specificity phosphatase 10 


2,27E-03 


1.99E-02 


-2,29 ^| 




MPST 


Mercaptopyruvate sulfurtransferase 


1 ,53E-02 


3,51 E-02 


-1,97 ■ 


Protein localization 


TMSB10 


Thymosin beta 10 


7,12E-03 


1,61 E-02 




Protein serine/threonine kinase activity 


lgFBp3 


Insulin-like growth factor binding protein 3 


3,38E-03 


4.03E-02 


-1,89 ■ 


Protein transport 


SFT2D2 


SFT2 domain containing 2 


1.50E-02 


3,26E-02 


-2,01 ■ 




SNX13 


Sorting nexin 13 


4,43E-02 


3.68E-02 


■ 1,95 




TMX4 


Thioredoxin-related transmembrane protein 4 


3.23E-02 


1.32E-02 


-2,50 


Proteolysis 


CTSK 


Cathepsin K 


1 ,67E-03 


4,85E-02 


-1,79 ■ 


Regulation of apoptosis 


NAPA 


N-ethylmaleimide-sensitive factor attachment protein, alpha 


1,41 E-02 


2.30E-02 


-2,21 | 


Transferase 


PRDM4 


PR domain containing 4 


3.99E-02 


4,83E-02 


-1,79 ■ 


Transmembrane region 


INSIG2 


Insulin induced gene 2 


1 ,72E-02 


1 ,78E-02 


-2,35 ^| 




KCTD1 


Potassium channel tetramerisation domain containing 1 


4,41 E-02 


3.66E-02 


-1,95 ■ 




0ST4 


Dolichyl-diphosphooligosaccharide-protein glycosyltransferase subunit 4 


3,27E-03 


3,05E-02 


■ 2,05 


Transport 


PFDN1 


Prefoldin subunit 1 


3,58E-02 


3.96E-02 


-1,90 ■ 




SLC44A1 


Solute carrier family 44, member 1 


4,85E-03 


4.03E-02 


■ 1,89 


Other 


BROX 


BR01 domain and CAAX motif containing 


3,26E-02 


3,32E-02 


-2,00 ■ 




CNOT10 


CCR4-NOT transcription complex, subunit 10 


3,63E-02 


1.46E-02 






C0PS4 


COP9 constitutive photomorphogenic homolog subunit 4 (Arabidopsis) 


1 ,93E-02 


4.12E-02 


-1,88 ■ 




DYNC1I2 


Dynein, cytoplasmic 1, intermediate chain 2 


2,39E-02 


3,24E-02 


■ 2,02 




GID8 


GID complex subunit 8 homolog (S. cerevisiae) 


1 ,63E-02 


1 ,54E-02 


-2,42 ^| 




KBTBD4 


Kelch repeat and BTB (POZ) domain containing 4 


4,37E-02 


4.47E-02 


-1,84 ■ 




KIAA2026 


KIAA2026 ortholog 


1 ,83E-02 


9.76E-03 


-2,66 




LOC 100037674 


Thymidylate synthetase 


2,60E-02 


3.98E-02 


-1,90 ■ 




LOC1 00847238 


Uncharacterized LOC1 00847238 


3,34E-02 


2,43E-02 


2 ,18 




LOC1 01 103408 


Uncharacterized LOC101 103408 


4,49E-02 


3,87E-02 


■ 1,92 




LOC1 01 105790 


Keratin-associated protein 15-1 -like 


3.97E-02 


4,62E-02 


-1,82 ■ 




LOC1 01 105958 


Non-histone chromosomal protein HMG-17-like 


9,15E-03 


3.40E-02 


-1,99 ■ 




LOC101108110 


Prothymosin alpha-like 


1 ,86E-02 


2,06E-03 


H 3 < 47 




LUU1 Ul \ uy /4/ 


Mamu class II histocompatibility antigen, DR alpha chain-like 


1.08E-03 


4,13E-03 


■ 3,11 




LOC101110664 


Cytochrome c oxidase subunit 5B, mitochondrial-like 


4,47E-03 


2,11 E-02 


-2,25 ^| 




LOC101 114533 


10 kDa heat shock protein, mitochondrial-like 


4,48E-02 


4,62E-02 


-1,82 ■ 




LOC101115083 


Glutathione S-transferase omega-1-like 


3,77E-03 


3.53E-02 


■ 1,97 




LOC101 117272 


Caspase-13-like 


1 ,28E-02 


2,94E-02 


-2,07 ■ 




LOC101119517 


Interferon-induced guanylate-binding protein 1 -like 


1.68E-02 


1,81 E-02 


-2,34 ^| 




LOC1 01 121538 


Cytochrome c oxidase subunit 6B1 -like 


1 ,35E-02 


8,38E-03 


2,74 




LOC101123112 


60S ribosomal protein L23a-like 


3,36E-02 


9,56E-03 


-2,67 




LOC523454 


Protein WWC3-like 


4,34E-02 


1.57E-02 


-2,41 ■ 




MY018A 


Myosin XVI MA 


2,30E-02 


4,81 E-02 


-1,79 ■ 




PAIP1 


Poly(A) binding protein interacting protein 1 


1 ,66E-02 


9,24E-03 


■ 2,69 




PPIA 


Peptidylprolyl isomerase A 


1 ,65E-02 


7.53E-03 


^| 2,80 




SECISBP2L 


SECIS binding protein 2-like 


1.36E-02 


2.15E-02 


-2,24 ^| 




SERBP1 


SERPINE1 mRNA binding protein 1 


3,82E-02 


1 ,94E-02 


■ 2,30 




SMU1 


Smu-1 suppressor of mec-8 and unc-52 homolog (C. elegans) 


4,77E-02 


1 ,38E-02 


2,48 




TSR2 


TSR2, 20S rRNA accumulation, homolog (S. cerevisiae) 


1.75E-02 


1 ,77E-02 


-2,35 ^| 



Figure 4 (See legend on next page.) 



Filali et al. BMC Genomics 2014, 15:59 
http://www.biomedcentral.eom/1 471 -21 64/1 5/59 



Page 15 of 17 



(See figure on previous page.) 

Figure 4 Relationship between gene expression profiles and PrP Sc deposition. Figure shows genes whose expression was statistically 
associated with PrP Sc deposition. P1 is the probability calculated using the simple regression model y, = p + bp-, + e„ and P2 is the probability 
calculated using the second regression model yy^n + Ti + bpy + ey (which includes a systematic effect associated with the 3 categories; Preclinical, 
Clinical and Healthy). The slope of regression describing the relationship between histopathological lesions and gene expression was calculated 
using the second model. Analyses were performed using R software (R project for Statistical Computing). 



expression and blood levels of CP will be required to de- 
termine its role in scrapie. 

We detected upregulation of gamma-butyrobetaine hy- 
droxylase (BBOX1), a gene involved in the response of 
the organism to oxidative stress. BBOX1 is the enzyme 
responsible for the biosynthesis of L-carnitine, a mol- 
ecule involved in fatty acid metabolism with important 
anti-free radical and antioxidant activities [66]. The acet- 
ylated form of L-carnitine (acetyl-L-carnitine) has been 
proposed as a potential treatment for neurodegenerative 
diseases, including AD [30]. Our results demonstrate up- 
regulation of BBOX1 in both clinical and preclinical 
stages of scrapie in the LRS, perhaps reflecting an in- 
crease in the levels of antioxidant molecules (L-carnitine 
and acetyl-L-carnitine) in response to the damage in- 
duced by the prion. 

The known functions of the genes that were up- and 
downregulated in scrapie suggest potential roles in the 
pathogenesis of prion diseases. However, we cannot rule 
out the possibility that these changes in expression are a 
secondary response of the organism to the disease. The 
findings of this transcriptomic analysis thus need to be 
complemented by further studies analyzing the cellular 
localization, protein levels and function of the genes of 
interest in order to establish their specific roles in prion 
diseases. 

Conclusions 

Our genome-wide expression analysis of the LRS of 
sheep with natural scrapie, at both preclinical and clin- 
ical stages of the disease, has identified new genes that 
may be involved in the pathogenesis of scrapie. We con- 
firmed the differential expression of 7 genes involved in 
prion or other neurodegenerative diseases. The down- 
regulation of genes involved in repairing misfolded pro- 
teins (PFDN2, PSMA7 and UCHL1) may contribute to 
prion formation. Furthermore, the altered expression of 
genes that promote angiogenesis in the lymph nodes 
(SERPINE1 and VEGFA) may facilitate prion dissemin- 
ation, while the observed decrease in CP expression may 
contribute to prion toxicity. Our results confirm the al- 
tered expression of BBOX1, a gene implicated in the 
systemic response of the organism to prion-induced oxi- 
dative stress, in scrapie infected animals. Although fur- 
ther studies are necessary to determine the specific roles 
of these genes in the pathogenesis of prion diseases, our 
findings represent an important first step towards 



identifying candidate genes that may serve as useful bio- 
markers for the development of early diagnostic tools 
and therapeutic strategies for scrapie. 

Availability of supporting data 

Microarray data are available in the ArrayExpress data- 
base: accession number E-MTAB-1346 (http://www.ebi. 
ac.uk/arrayexpress/experiments/E-MTAB-1346/). 

Abbreviations 

ARQ: Alanine arginine glutamine; BB0X1: Gamma-butyrobetaine hydroxylase; 
BSE: Bovine spongiform encephalopathy; CNS: Central nervous system; 
CP: Ceruloplasmin; DAB: 3,3'-diaminobenzidine; FDCs: Follicular dendritic 
cells; G6PDH: Glucose-6-phosphate dehydrogenase; GAPDH: Glyceraldehyde- 
3-phosphate dehydrogenase; GO: Gene ontology; IHC: Immunohistochemistry; 
LRS: Lymphoreticular system; MLN: Mesenteric lymph node; PFDN2: Prefoldin 
subunit 2; PRNP: Prion protein; PrP c : Cellular prion protein; PrP Sc : Scrapie prion 
protein; PSMA7: Proteasome subunit alpha type-7; RIN: RNA integrity number; 
RPL32: Ribosomal protein I32; SERPINE1: Plasminogen activator inhibitor-!; 
TBMs: Tingible body macrophages; TSEs: Transmissible spongiform 
encephalopathies; UCHL1: Carboxy-terminal hydrolase L1; VEGFA: Vascular 
endothelial growth factor; VRQ: Valine arginine glutamine. 

Competing interests 

The authors declare that they have no competing interests. 
Authors' contributions 

HF performed the experiments and drafted the manuscript. 1MB helped 
design the study, participated in the molecular genetic studies and the 
sequence alignment study, and drafted the manuscript. FH participated in 
the microarray analysis. LV performed the association analysis. CH 
participated in the pathological characterization of the animals. DM 
performed the qRT-PCR experiments. MM helped design the study and 
drafted the manuscript. AB helped design the microarray and sequence 
alignment studies. JJB helped design and coordinate the study and to draft 
the manuscript. RB conceived the study, participated in its design and 
coordination, and helped to draft the manuscript. All authors read and 
approved the final version of the manuscript. 

Acknowledgements 

The authors wish to acknowledge the excellent technical support of Belen 
Mann, Silvia Ruiz, Nuria Segovia, Silvia Castejon, and Sonia Gomez. This work was 
performed as part of the AGL2008-0256 project, financed by MICINN-FEDER. 
Hicham Filali was supported by doctoral grants from the Ministerio de Asuntos 
Exteriores y de Cooperacion (MAEC/AECID). The arrays and sequenced libraries 
were made available from projects at CVI supported by grants from the Dutch 
Ministry of Agriculture, Nature Management, and Fisheries (LNV). 

Author details 

Centra de Investigacion en Encefalopatias y Enfermedades Transmisibles 
Emergentes, Facultad de Veterinaria, Universidad de Zaragoza, Zaragoza, 
Spain. 2 Laboratorio de Genetica Bioquimica (LAGENBIO), Facultad de 
Veterinaria, Universidad de Zaragoza, Zaragoza, Spain. 3 Central Veterinary 
Institute part of Wageningen UR (CVI), Lelystad, The Netherlands. 4 Unidad de 
Genetica Cuantitativa y Mejora Animal, Facultad de Veterinaria, Universidad 
de Zaragoza, Zaragoza, Spain. 

Received: 24 January 2013 Accepted: 17 January 2014 
Published: 23 January 2014 



Filali et al. BMC Genomics 2014, 15:59 
http://www.biomedcentral.eom/1 471 -21 64/1 5/59 



Page 16 of 17 



References 

1. Aguzzi A, Heikenwalder M: Pathogenesis of prion diseases: current status 

and future outlook. Not Rev Microbiol 2006, 4:765-775. 

2. Prusiner SB: Prions. Proc Natl Acad Sci USA 1 998, 95:1 3363-1 3383. 

3. Hunter N: Scrapie and experimental BSE in sheep. B Med Bull 2003, 
66:171-183. 

4. Jeffrey M, Gonzalez L: Classical sheep transmissible spongiform 
encephalopathies: pathogenesis, pathological phenotypes and clinical 
disease. Neuropathol App Neurobiol 2007, 33:373-394. 

5. Andreoletti O, Berthon P, Marc D, Sarradin P, Grosclaude J, Van Keulen L, 
Schelcher F, Elsen J, Lantier F: Early accumulation of PrPsc in gut- 
associated lymphoid and nervous tissues of susceptible sheep from a 
Romanov flock with natural scrapie. J Gen Virol 2000, 81:31 15-3126. 

6. Van Keulen LJ, Vromans ME, Van Zijderveld FG: Early and late pathogenesis 
of natural scrapie infection in sheep. APMIS 2002, 1 10:23-32. 

7. Van Keulen LJ, Bossers A, Van Zijderveld F: TSE pathogenesis in cattle and 
sheep. Vet Res 2008, 39:24. 

8. Brown KL, Stewart K, Ritchie DL, Mabbott NA, Williams A, Fraser H, Morrison 
Wl, Bruce ME: Scrapie replication in lymphoid tissues depends on prion 
protein-expressing follicular dendritic cells. Nat Med 1999, 5:1308-1312. 

9. McGovern G, Jeffrey M: Scrapie-specific pathology of sheep lymphoid 
tissues. PloS One 2007, 2:e1304. 

10. Filali H, Martin-Burriel I, Harders F, Varona L, Lyahyai J, Zaragoza P, Pumarola 
M, Badiola JJ, Bossers A, Bolea R: Gene expression profiling and association 
with prion-related lesions in the medulla oblongata of symptomatic 
natural scrapie animals. PloS One 201 1, 6:e19909. 

1 1 . Filali H, Martin-Burriel I, Harders F, Varona L, Serrano C, Acin C, Badiola JJ, 
Bossers A, Bolea R: Medulla oblongata transcriptome changes during 
presymptomatic natural scrapie and their association with prion-related 
lesions. BMC Genomics 2012, 13:399. 

12. Skinner PJ, Abbassi H, Chesebro B, Race RE, Reilly C, Haase AT: Gene 
expression alterations in brains of mice infected with three strains of 
scrapie. BMC Genomics 2006, 7:1 14. 

13. Brown AR, Webb J, Rebus S, Williams A, Fazakerley JK: Identification of up- 
regulated genes by array analysis in sera pie-infected mouse brains. 
Neuropathol App Neurobiol 2004, 30:555-567. 

14. Xiang W, Windl O, Wunsch G, Dugas M, Kohlmann A, Dierkes N, Westner IM, 
Kretzschmar HA: Identification of differentially expressed genes in 

sera pie-infected mouse brains by using global gene expression technology. 
J Virol 2004, 78:11051-11060. 

1 5. Guillerme-Bosselut F, Forestier L, Jayat-Vignoles C, Vilotte JL, Popa I, Portoukalian 
J, Le Dur A, Laude H, Julien R, Gallet PF: Glycosylation-related gene expression 
profiling in the brain and spleen of sera pie-affected mouse. Glycobiology 
2009,19:879-889. 

16. Huzarewich RL, Medina S, Robertson C, Parchaliuk D, Booth SA: 
Transcriptional modulation in a leukocyte-depleted splenic cell 
population during prion disease. J Toxicol Environ Health Part A 201 1, 
74:1504-1520. 

17. Skretting G, Espenes A, Ulvund MJ, Olsaker I: cDNA representational 
difference analysis of ileal Peyer's patches in lambs after oral inoculation 
with scrapie. Biochem Biophys Res Commun 2004, 316:272-279. 

18. Vargas F, Lujan L, Bolea R, Monleon E, Martin-Burriel I, Fernandez A, De Bias 
I, Badiola JJ: Detection and clinical evolution of scrapie in sheep by 3rd 
eyelid biopsy. J Vet Intern Med 2006, 20:187-193. 

19. Monleon E, Garza MC, Sarasa R, Alvarez-Rodriguez J, Bolea R, Monzon M, 
Vargas MA, Badiola JJ, Acin C: An assessment of the efficiency of PrPsc 
detection in rectal mucosa and third-eyelid biopsies from animals 
infected with scrapie. Vet Microbiol 201 1, 147:237-243. 

20. Monleon E, Monzon M, Hortells P, Vargas A, Acin C, Badiola JJ: Detection of 
PrPsc on lymphoid tissues from naturally affected scrapie animals: 
comparison of three visualization systems. J Histochem Cytochem 2004, 
52:145-151. 

21. Caraux G, Pinloche S: PermutMatrix: a graphical environment to arrange 
gene expression profiles in optimal linear order. Bioinformatics 2005, 
21:1280-1281. 

22. Dennis G Jr, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, Lempicki RA: 
DAVID: database for annotation, visualization, and integrated discovery. 
Genome Biol 2003, 4:P3. 

23. Da Huang W, Sherman BT, Lempicki RA: Systematic and integrative 
analysis of large gene lists using DAVID bioinformatics resources. 
Nat Protoc 2009, 4:44-57. 



24. Del Bo R, Ghezzi S, Scarpini E, Bresolin N, Comi GP: VEGF genetic variability 
is associated with increased risk of developing Alzheimer's disease. 

J Neurol Sci 2009, 283:66-68. 

25. Fang Y, Zhang L, Zeng Z, Lian Y, Jia Y, Zhu H, Xu Y: Promoter 
polymorphisms of SERPINE1 are associated with the antidepressant 
response to depression in Alzheimer's disease. Neurosci Lett 2012, 
516:217-220. 

26. Mays CE, Ryou C: Plasminogen stimulates propagation of protease- 
resistant prion protein in vitro. FASEB J 2010, 24:5102-51 12. 

27. Saigoh K, Wang YL, Suh JG, Yamanishi T, Sakai Y, Kiyosawa H, Harada T, 
Ichihara N, Wakana S, Kikuchi T, Wada K: Intragenic deletion in the gene 
encoding ubiquitin carboxy-terminal hydrolase in gad mice. Nat Genet 
1999, 23:47-51. 

28. Deriziotis P, Tabrizi SJ: Prions and the proteasome. Biochim Biophys Acta 
2008, 1782:713-722. 

29. Morita H, Ikeda S, Yamamoto K, Morita S, Yoshida K, Nomoto S, Kato M, 
Yanagisawa N: Hereditary ceruloplasmin deficiency with hemosiderosis: a 
clinicopathological study of a Japanese family. Ann Neurol 1995, 
37:646-656. 

30. Aliev G, Li Y, Palacios HH, Obrenovich ME: Oxidative stress induced 
mitochondrial DNA deletion as a hallmark for the drug development in 
the context of the cerebrovascular diseases. Recent Pat Cardiovasc Drug 
Discov 2011,6:222-241. 

31. Lyahyai J, Serrano C, Ranera B, Badiola JJ, Zaragoza P, Martin-Burriel I: Effect 
of scrapie on the stability of housekeeping genes. Anim Biotechnol 2010, 
21:1-13. 

32. Garcia-Crespo D, Juste RA, Hurtado A: Selection of ovine housekeeping 
genes for normalisation by real-time RT-PCR; analysis of PrP gene 
expression and genetic susceptibility to scrapie. BMC Vet Res 2005, 1:3. 

33. Gossner A, Roupaka S, Foster J, Hunter N, Hopkins J: Transcriptional 
profiling of peripheral lymphoid tissue reveals genes and networks 
linked to SSBP/1 scrapie pathology in sheep. Vet Microbiol 201 1, 
153:218-228. 

34. Kim HO, Snyder GP, Blazey TM, Race RE, Chesebro B, Skinner PJ: Prion 
disease induced alterations in gene expression in spleen and brain prior 
to clinical symptoms. Adv Appl Bioinformatics Chem: AABC 2008, 1:29-50. 

35. Julius C, Hutter G, Wagner U, Seeger H, Kana V, Kranich J, Klohn PC, 
Weissmann C, Miele G, Aguzzi A: Transcriptional stability of cultured cells 
upon prion infection (vol 375, pg 1222, 2008). J Mol Biol 2009, 
388:207-207. 

36. Stanton JB, Knowles DP, Call DR, Mathison BA, Baszler TV: Limited 
transcriptional response of ovine microglia to prion accumulation. 
Biochem Biophys Res Commun 2009, 386:345-350. 

37. Fasano C, Campana V, Griffiths B, Kelly G, Schiavo G, Zurzolo C: Gene 
expression profile of quinacrine-cured prion-infected mouse neuronal 
cells. J Neurochem 2008, 105:239-250. 

38. Booth S, Bowman C, Baumgartner R, Sorensen G, Robertson C, Coulthart M, 
Phillipson C, Somorjai RL: Identification of central nervous system genes 
involved in the host response to the scrapie agent during preclinical 
and clinical infection. J Gen Virol 2004, 85:3459-3471. 

39. Vidal E, Acin C, Foradada L, Monzon M, Marquez M, Monleon E, Pumarola 
M, Badiola JJ, Bolea R: Immunohistochemical characterisation of classical 
scrapie neuropathology in sheep. J Comp Pathol 2009, 141:135-146. 

40. Costa C, Tortosa R, Domenech A, Vidal E, Pumarola M, Bassols A: Mapping 
of aggrecan, hyaluronic acid, heparan sulphate proteoglycans and 
aquaporin 4 in the central nervous system of the mouse. J Chem 
Neuroanat 2007, 33:111-123. 

41 . Loubet D, Dakowski C, Pietri M, Pradines E, Bernard S, Callebert J, Ardila-Osorio 
H, Mouillet-Richard S, Launay JM, Kellermann O, Schneider B: Neuritogenesis: 
the prion protein controls betal integrin signaling activity. FASEB J 2012, 
26:678-690. 

42. Kolsch H, Lehmann DJ, Ibrahim-Verbaas CA, Combarros O, Van Duijn CM, 
Hammond N, Belbin O, Cortina-Borja M, Lehmann MG, Aulchenko YS, et al: 
Interaction of insulin and PPAR-alpha genes in Alzheimer's disease: the 
Epistasis Project. J Neural Transm 2012, 1 19:473-479. 

43. Mandrekar-Colucci S, Karlo JC, Landreth GE: Mechanisms underlying the 
rapid peroxisome proliferator-activated receptor-gamma-mediated 
amyloid clearance and reversal of cognitive deficits in a murine model 
of Alzheimer's disease. J Neurosci 201 2, 32:1 01 1 7-1 01 28. 

44. Dandoy-Dron F, Guillo F, Benboudjema L, Deslys JP, Lasmezas C, Dormont 
D, Tovey MG, Dron M: Gene expression in scrapie: cloning of a new 



Filali et al. BMC Genomics 201 4, 1 5:59 Page 1 7 of 1 7 

http://www.biomedcentral.eom/1 471 -21 64/1 5/59 



45. 



47. 



49. 



50. 



52. 



53. 



54. 



55. 



56. 



57. 



58. 



59. 



60. 



62. 



63. 



64. 



65. 



66. 



scrapie-responsive gene and the identification of increased levels of 
seven other mRNA transcripts. J Biol Chem 1998, 273:7691-7697. 
Baker CA, Manuelidis L: Unique inflammatory RNA profiles of microglia in 
Creutzfeldt-Jakob disease. Proc Natl Acad Sci USA 2003, 100:675-679. 
Nakamura Y, Takeda M, Yoshimi K, Hattori H, Hariguchi S, Kitajima S, 
Hashimoto S, Nishimura T: Involvement of clathrin light chains in the 
pathology of Alzheimer's disease. Acta Neuropathol 1994, 87:23-31. 
Legendre-Guillemin V, Metzler M, Lemaire JF, Philie J, Gan L, Hayden MR, 
McPherson PS: Huntingtin interacting protein 1 (HIP1) regulates clathrin 
assembly through direct binding to the regulatory region of the clathrin 
light chain. J Biol Chem 2005, 280:6101-6108. 

Lambert JC, Amouyel P: Genetics of Alzheimer's disease: new evidences 
for an old hypothesis? Curr Opin Genet Dev 201 1, 21:295-301. 
Cook C, Petrucelli L: A critical evaluation of the ubiquitin-proteasome 
system in Parkinson's disease. Biochim Biophys Acta 2009, 1792:664-675. 
Goldberg AL: Protein degradation and protection against misfolded or 
damaged proteins. Nature 2003, 426:895-899. 

Wilkinson KD, Lee KM, Deshpande S, Duerksen-Hughes P, Boss JM, Pohl J: 
The neuron-specific protein PGP 9.5 is a ubiquitin carboxyl-terminal 
hydrolase. Science 1989, 246:670-673. 

Choi J, Levey Al, Weintraub ST, Rees HD, Gearing M, Chin LS, Li L: Oxidative 
modifications and down-regulation of ubiquitin carboxyl-terminal 
hydrolase L1 associated with idiopathic Parkinson's and Alzheimer's 
diseases. J Biol Chem 2004, 279:13256-13264. 

Naze P, Vuillaume I, Destee A, Pasquier F, Sablonniere B: Mutation analysis 
and association studies of the ubiquitin carboxy-terminal hydrolase L1 
gene in Huntington's disease. Neurosci Lett 2002, 328:1-4. 
Moore DJ, Dawson VL, Dawson TM: Role for the ubiquitin-proteasome 
system in Parkinson's disease and other neurodegenerative brain 
amyloidoses. Neuromol Med 2003, 4:95-108. 

Dawson TM, Dawson VL: Molecular pathways of neurodegeneration in 
Parkinson's disease. Science 2003, 302:819-822. 

Deriziotis P, Andre R, Smith DM, Goold R, Kinghorn KJ, Kristiansen M, Nathan 
JA, Rosenzweig R, Krutauz D, Glickman MH, et al: Misfolded PrP impairs the 
UPS by interaction with the 20S proteasome and inhibition of substrate 
entry. EMBO J 201 1, 30:3065-3077. 

Ma J, Lindquist S: Conversion of PrP to a self-perpetuating PrPSc-like 
conformation in the cytosol. Science 2002, 298:1785-1788. 
Kang SC, Brown DR, Whiteman M, Li R, Pan T, Perry G, Wisniewski T, Sy MS, 
Wong BS: Prion protein is ubiquitinated after developing protease 
resistance in the brains of scrapie-infected mice. J Pathol 2004, 203:603-608. 
Vainberg IE, Lewis SA, Rommelaere H, Ampe C, Vandekerckhove J, Klein HL, 
Cowan NJ: Prefoldin, a chaperone that delivers unfolded proteins to 
cytosolic chaperonin. Cell 1998, 93:863-873. 
Broer L, Ikram MA, Schuur M, DeStefano AL, Bis JC, Liu F, Rivadeneira F, 
Uitterlinden AG, Beiser AS, Longstreth WT, et al: Association of HSP70 and its 
co-chaperones with Alzheimer's disease. J Alzheimers Dis 201 1, 25:93-102. 
Bajou K, Maillard C, Jost M, Lijnen RH, Gils A, Declerck P, Carmeliet P, Foidart 
JM, Noel A: Host-derived plasminogen activator inhibitor-1 (PAI-1) 
concentration is critical for in vivo tumoral angiogenesis and growth. 
Oncogene 2004, 23:6986-6990. 

Hirakawa S, Kodama S, Kunstfeld R, Kajiya K, Brown LF, Detmar M: VEGF-A 
induces tumor and sentinel lymph node lymphangiogenesis and 
promotes lymphatic metastasis. J Exp Med 2005, 201:1089-1099. 
Kincaid AE, Hudson KF, Richey MW, Bartz JC: Rapid transepithelial 
transport of prions following inhalation. J Virol 2012, 86:12731-12740. 
Waggoner DJ, Bartnikas TB, Gitlin JD: The role of copper in 
neurodegenerative disease. Neurobiol Dis 1999, 6:221-230. 
Dexter DT, Jenner P, Schapira AH, Marsden CD: Alterations in levels of iron, 
ferritin, and other trace metals in neurodegenerative diseases affecting 
the basal ganglia: the Royal Kings and Queens Parkinson's Disease 
Research Group. Ann Neurol 1992, 32(Suppl):S94-S100. 
Gulcin I: Antioxidant and antiradical activities of L-carnitine. Life Sci 2006, 
78:803-811. 



doi:1 0.1 1 86/1 471 -21 64-1 5-59 

Cite this article as: Filali et al.: Gene expression profiling of mesenteric 
lymph nodes from sheep with natural scrapie. BMC Genomics 201 4 15:59. 



Submit your next manuscript to BioMed Central 
and take full advantage of: 

• Convenient online submission 

• Thorough peer review 

• No space constraints or color figure charges 

• Immediate publication on acceptance 

• Inclusion in PubMed, CAS, Scopus and Google Scholar 

• Research which is freely available for redistribution 



Submit your manuscript at 
www.biomedcentral.com/submit 



o 



BioMed Central 



